Thermoelectric effects in superconducting proximity structures 
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Attaching a superconductor in good contact with a normal metal makes rise to a proximity effect where the 
superconducting correlations leak into the normal metal. An additional contact close to the first one makes 
it possible to carry a supercurrent through the metal. Forcing this supercurrent flow along with an additional 
quasiparticle current from one or many normal-metal reservoirs makes rise to many interesting effects. The 
supercurrent can be used to tune the local energy distribution function of the electrons. This mechanism also 
leads to finite thermoelectric effects even in the presence of electron-hole symmetry. Here we review these 
effects and discuss to which extent the existing observations of thermoelectric effects in metallic samples can 
be explained through the use of the dirty-limit quasiclassical theory. 



I. INTRODUCTION 

Applying a bias voltage or a temperature gradient across a 
conductor makes rise to charge and energy currents. The lin- 
ear response between the biases and currents is described via 
the thermoelectric matrix, whose diagonal parts are the charge 
and thermal conductances, and the off-diagonal parts are often 
referred to as the thermoelectric coefficients. In typical met- 
als, the latter arise due to the asymmetry between positive- and 
negative-energy excitations with respect to the Fermi energy, 
i.e., electrons and holes. Such asymmetry in metals is very 
small, making the typical thermoelectric effects at sub-Kelvin 
temperatures hard to measure accurately. 

Placing a superconductor in good contact to a normal-metal 
conductor makes rise to finite pair correlations also inside the 
latter, even when the pair potential inside it vanishes. This su- 
perconducting proximity effect has an energy-dependent pen- 
etration depth; at typical measurement temperatures of the or- 
der of 100 mK it extends up to the micrometer range. The 
proximity effect modifies the thermoelectric response of the 
normal conductor. Most importantly, it makes rise to ther- 
moelectric effects which are orders of magnitude larger than 
in the absence of superconductivity. The proximity-induced 
modifications are discussed in this paper by employing the 
quasiclassical theory in the diffusive limi t i 2 ' 3 . In this theory, 
we assume that all the relevant length scales of the problem 
exceed especially the Fermi wavelength (quasiclassical ap- 
proximation) and the mean free path (diffusive limit). Exam- 
ples of such relevant length scales are the structure size and the 
superconducting coherence length. A further property of the 
quasiclassical theory, especially important for thermoelectric 
effects, is that it assumes electron-hole symmetry. Because of 
this, in the normal state it predicts vanishing thermoelectric 
coefficients. 

The proximity modification of the thermoelectric matrix 
is conveniently described in Andreev interferometers (see 
Fig. [TJ, where there are two superconducting contacts to the 
normal metal. In this structure, the phase difference between 
the two contacts affects the proximity modifications, and its 
presence is an important requirement for finite thermoelectric 
effects, at least within the quasiclassical theory. This type of 
a dependence of the electric conductance on the phase has for 
example been suggested for use in quantum measurements of 



flux qubitsi 

This paper is organized as follows. In Sec. [II] we briefly in- 
troduce the thermoelectric effects and their relations in normal 
metals and after that detail the quasiclassical equations for the 
diffusive limit. We also discuss briefly effects left out in the 
present paper. Section [HI] introduces to the properties of the 
spectral supercurrent and two aspects of nonequilibrium su- 
percurrent: first, we briefly mention the nonequilibrium con- 
trol of the supercurrent and then concentrate more on how the 
energy distribution function of electrons is controlled with the 
supercurrent. The latter is a precursor to the thermoelectric 
effects. These are described in Sec. lIVI which first introduces 
to the special symmetries of the thermoelectric matrix, and 
then details the behavior of its components as a function of 
the phase difference and temperature. In Sec.|V] we mention 
the effects relevant when considering the dependence on the 
magnetic flux, and finally in Sec. [VT]we conclude and point 
out the open questions related with understanding the mea- 
surements of the thermoelectric effects. 



II. THERMOELECTRIC TRANSPORT IN PROXIMITY 
STRUCTURES 

In this paper, we show that the superconducting proximity 
effect is able to generate large thermoelectric effects, which 
can be described without employing the electron-hole asym- 
metry. An important factor in the theory is the presence of 
supercurrent, which then needs to be taken into account in 
the description of the currents. Moreover, to describe trans- 
port between normal metals in the presence of supercurrent, 
we need to have multiple terminals connected to the structure. 
As the biases mostly deal with the quasiparticle current, we 
define the thermoelectric matrix in a multi-terminal structure 
according to 
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where I' c and I' s eq are the total charge current and the equi- 
librium supercurrent flowing to terminal i, and I' E is the en- 
ergy current (supercurrent as such carries no energy current). 
Moreover, AVj is the bias voltage and ATj is the temperature 
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difference from some average temperature T, both present in 
terminal j. 



A. Transport in normal-metal structures 

Thermoelectricity in normal-metal wires can be practically 
described especially in the diffusive limit (structure size L, 
elastic mean free path £ e \ and Fermi wavelength Xf satisfying 
the relation L > { e i > Xf). In this limit, the charge and heat 
currents flowing in the wire are given by 

I c = —eA J dED{E)v(E)d x f(x;E) (2a) 

Iq = -A [ dE(E - n)D{E)v{E)d x f{x-E). (2b) 

Here D(E) is the diffusion constant, v(E) is the density of 
states, A is the cross-sectional area of the wire and x is the 
coordinate parallel to the wire. The heat current can be simply 
related to the energy current via Iq=I% — fll c - For linear 
response, Iq = Ie — the second term is responsible for Joule 
heating — and we can expand the electron energy distribution 
function d x f = (drf)d x T + (<9jj/)<9 x /x, with / fy, the Fermi 
function. Furthermore, assuming some characteristic length L 
and taking d x T = AT /L and d x ji = eAV /L allows us to relate 
the results to Eq. (Q~|). 

The energy-dependent changes in the density of states or 
the diffusion constant typically take place at large energy 
scales of the order of the Fermi energy Ep . We can thus 
expand them as D(E) i=s D + cq{E — Ef)/Ep and v(E) « 
Vf + c n{E — Ef)/Ef. To linear order in cjy and c^, we then 
ge£ 

L\\ = G = e 2 VpDA/L, Drude conductance (3a) 

L22=LqGT 2 , Wiedemann-Franz law (3b) 

L 12 = eLoG'T 2 , Mott law (3c) 

L21 = L\2, Onsager-Kelvin relation. (3d) 

Here L = n 2 k 2 B /(3e 2 ) « 2.45 x 10" 8 Wi2K~ 2 is the Lorenz 
number, and the electron-hole asymmetry is described by the 
factor G' — e 2 (coVF + Dcn)A/(LEf). These relations show 
that the thermoelectric effects in normal metals are of the or- 
der of kgT jEp. 

The Onsager-Kelvin relation between the two thermoelec- 
tric coefficients is an example of a more general relation 6,7 ' 8 
between different linear-response coefficients. According to 
this relation, the elements of the thermoelectric matrix in 
Eq. ((TJ should satisfy 

under the reversal of the magnetic field B. Here a, /3 € {1,2}. 
This relation results essentially only from the assumption of 
time-reversal symmetry. In Sec. lIV Bl we show how this equa- 
tion can be derived for the energy-dependent response coeffi- 
cients within the quasiclassical theory. 



The presence of superconductivity modifies the above 
laws in many different ways 9 . For example, the Andreev 
reflection 10 breaks the Wiedemann-Franz law, and the Mott 
law is broken in asymmetric structures 11 . The effects re- 
lated to the superconducting density of states or to charge 
imbalance make modifications to the thermoelectric effects 
at interface s 12 ' 13 and for the nonlinear response 14 . The 
main modification at linear response due to the proximity ef- 
fect is the appearance of thermoelectric effects even without 
electron-hole asymmetry 1 5 - 16 i 1 7 - 1 8 ' 1 9 ' 2Q . The latter effect is 
at the low temperatures where superconductivity can be ob- 
served much stronger than that expected from the electron- 
hole asymmetry. Therefore, we concentrate on an electron- 
hole symmetric theory in the remainder of this paper. We em- 
ploy the quasiclassical theory that provides a fair description 
of inhomogeneous superconductivity both in equilibrium and 
nonequilibrium systems. Moreover, for simplicity and also 
dictated by many of the experiments, we concentrate on the 
diffusive limit. 

B. Usadel equations for proximity structures 

Heterostructures composed of diffusive normal-metal or su- 
perconducting wires in and out of equilibrium can be de- 
scribed through the use of Usadel equations 1 for the Keldysh 
Green's functions G. These equations are reviewed in many 
references — we cite here only a few of those 218 apply- 
ing similar parametrization as here. Written in the Nambu 
<g> Keldysh space, Usadel equation is a nonlinear differential 
equation for a 4 x 4 matrix, 

0[Y,G[V,G]] = [-iE + A + t,G] . (5) 

Here D is the diffusion constant, E is the energy calcu- 
lated from the Fermi energy, A denotes the superconducting 
order parameter and E the self-energy for inelastic scatter- 
ing (mainly the part of electron-electron interaction not de- 
scribed by A and electron-phonon scattering), for spin-flip 
or spin-orbit scattering. In the presence of a magnetic field, 
V = V — /eAT3 is the gauge-invariant derivative including the 
vector potential A. In addition to Eq. ©, G satisfies the nor- 
malization G 2 = 1, where 1 is the identity matrix. 

In the diffusive limit, we implicitly assume that all the 
length scales of the problem, including the superconducting 
coherence length and the mean free paths for other types of 
scattering than elastic, are much longer than the elastic mean 
free path. An example of such other types of scattering is the 
spin-flip scattering, described in the Born approximation by 
the self-energy 2 -^ 

Ev/' = 7T— %Gf 3 , (6) 
2% 

where tw is the spin-flip scattering time. This term is included 
in the following analytic expressions, but omitted from the 
numerics. 

In the Keldysh space, Green's function has the form 

^_(G R G K \ 
G -{o &)> 
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where G R I A I K denote the Retarded/ Advanced/Keldysh func- 
tions. The latter are 2 x 2 matrices in the Nambu particle- 
hole space. 83 Products of this type of matrices yield similar 
matrices, without mixing the Keldysh parts into the diago- 
nal. Therefore, also the Usadel equation (0 has a similar ma- 
trix structure. Employing the normalization and the symmetry 
G A = — tiG R Ti, we may parametrize 



and 



G R = cosh(0)T3 + sinh(0)(cos(^)/T2 + sin(#)z'ti) 



G K = G R (f L +f T T 3 ) - (f L + f T r 3 )d A 



Here and % are complex scalar parameters, roughly describ- 
ing the magnitude and phase of the pair amplitude, respec- 
tively. In the Keldysh part, the additional parameters fi and 
fr are the longitudinal and transverse parts of the electron dis- 
tribution function. It can be shown 21 that this parametrization 
spans all the possible solutions of the Keldysh-Usadel equa- 
tions in non-magnetic systems. 
Usadel equations for and % are 

DV 2 9 = - 2i(E + 2'r in ) sinh 9 + (— + ^] sinh(20) 

W 2D J 



-2?'|A|cos(0 -%)cosh(0), 



V- (-v 5 sinh 2 0) = - 2z|A| sin(0 -£)sinh(0). 



v s =D(yX-2eA/h). 



(7a) 
(7b) 



Here we assume that the superconducting order parameter is 
of the form A = |A|e'^. Note that in a proximity structure, the 
superfluid velocity v$ is position dependent. We include the 
effect of weak inelastic scattering through a constant imagi- 
nary part Tin of the energy 22 . In the numerics, this is set to a 
small but finite positive value in order to preserve the analytic 
structure of the Green's functions. 

The kinetic equations for the distribution functions read 



DV-r T f=(V-j s )f L + 2\A\<%f T , 
fY/ = 3rVfr + ^ h + Ml , 

DV-fV = 0, 

t L f=3) L Vf L -3Vf T +j s f T , 



(8a) 
(8b) 



where the kinetic coefficients are 





i(l + |cosh0| 2 -|sinh0 


2 cosh(2Im[^])), 


(9a) 


&r = 


^(l + |cosh0| 2 + |sinh0 


2 cosh(2Im[ Z ])), 


(9b) 




i|sinh0| 2 sinh(2ImM)), 




(9c) 


h = 


lm[-sinh 2 (0)^], 




(9d) 




lm[-cos(0 -%)sinh(0)] 




(9e) 



Inside a superconductor where the pair interaction parameter 
A^O, the superconducting pair potential is obtained via 

A = - dE\ (e ix sinh d+e i% * sinh 9*)f L 

4 J ,J (10) 

- (e ix sinh - e ix ' sinh 9*)f T ] 

Solving Eqs. (O, © and ( TTOb we get the observables, for ex- 
ample the charge and energy current densities given by 

J c = -?-J~JEt T f, J E = ^J~JEEt L f. (11) 

In most of the text below, we assume that the superconductors 
are bulky reservoirs, such that the self-consistency equation 
([Tol l can be ignored. We rather concentrate on the phenomena 
taking place in normal-metal wires close to the superconduc- 
tors. In those wires, we assume A = 0, and thereby also A = 0. 
This simplifies the resulting equations. 



C. Interfaces and terminals 

Usadel equation holds within the wires where changes in 
the parameters take place slowly compared to the mean free 
path. At interfaces, it has to be supplemented by boundary 
conditions. Initially, these were derived for a general quasi- 
classical Green's function by Zaitsev 23 . For the diffusive case, 
the general boundary conditions were solved by Nazarov 24 . 
They read 



h 



2e 2 y T„[G Ll G R ] 
7Ch^4-T n ({G L ,G R }-2y (12) 
OiAtGiVGi ■ n, 



where a and A are the normal-state conductivity and the cross 
section of the wires next to the interface, subscript L/R denote 
left/right from the interface, and n is the unit vector perpen- 
dicular to the interface, pointing to the right. The interface 
is characterized by the set {t„} of transmission eigenvalues. 
Note that the resulting expression for / is linear in the electron 
distribution functions /, due to the Keldysh block structure of 
the Green's functions. In what follows, we assume that the 
normal-state conductance Gj = 2e 2 £„ X n /h for each interface 
is large, such that their effect can be neglected. However, the 
arguments on the general symmetries of the thermoelectric co- 
efficients are independent of this assumption. 

In addition to having the correct boundary conditions for in- 
terfaces, one needs also to describe the behavior of the Green's 
function G inside different types of terminals. A typical as- 
sumption is that the Green's functions obtain their bulk values 
very close to the interface between a wire and a terminal. Es- 
sentially this means that the specific resistance (both charge 
and thermal) of the terminals should be much smaller than 
that of the mesoscopic region under study. Experimentally 
this is realized by making the cross section of especially the 
normal-metal terminals much larger than that of the wires. 

Inside superconductors for energies E < |A| all quantities 
except fi relax to their bulk values wi thin dista nces compa- 
rable to the coherence length <^o = \JhD/ (2A). Moreover, 
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for the L-mode, the Andreev reflection boundary condition 
Yif = applies at these energies. These details of NS in- 
terfaces one can usually describe by increasing the effective 
length 25 of the normal-metal wires in contact to superconduc- 
tors by an amount comparable to <jjo- However, at energies 
E > |A|, nonequilibrium in fj and fi may persist to greater 
distances. This charge and energy imbalance is limited by in- 
elastic relaxation processes, and for the charge mode in the 
diffusive limit, by the decoherence induced by a flowing su- 
percurrent or spin-flip scattering. (See for example Refs . 22 > 26 .) 

For temperatures or voltages of the order or larger than A, 
we hence have to pay some attention to a proper treatment 
of superconductors, especially superconducting loops (for an 
example, see Fig. [TJ with length Le, cross section Al and 
normal-state conductivity Ge- Assume such a loop is con- 
nected to a normal-metal wire with length L w , cross section 
A w and normal-state conductivity a w . When compared to the 
superconductor, the latter is described by an effective length 
L' w = L w OlAl/ '(<J w A w ) to account for the differences in the 
specific resistance. Furthermore, assume an energy relaxation 
length Le inside the superconductor. We then have three prac- 
tically important limits: a) Le <§; Ll,L' w , b) L\ v <C Le <C Le 
and c) L' w <C Le <C Ll- In the first case, the relaxation in the 
superconductors is fast, and we may assume that fjSE > |A|) 
and fr{E > |A|) acquire their bulk values immediately at the 
superconducting interface. In the case b), the normal-state 
resistance of the loop is much higher than that of the normal- 
metal wires, so that the proper boundary condition is the van- 
ishing of quasiparticle current to the superconductors. For the 
case c), we again get a vanishing of the quasiparticle charge 
current, but the energy current will depend on the details of 
inelastic relaxation in the superconductor. 

We see no way to formulate exact mathematical boundary 
conditions for the limits b) and c) above — they in principle re- 
quire the solution of the Usadel equation inside the supercon- 
ductor. One attempt to approximate the case b) in a way con- 
sistent with the Onsager symmetry is described in Sec. |IVB"Tl 
It captures most of the essential physics of this problem, i.e., 
taking into account the finite charge and thermal resistance of 
the loop at high temperatures. 



D. Effects left out 

There are two more practically important self-energies that 
were not included in the above description: those related to 
electron-electron and electron-phonon interactions, E e -e and 
£e-ph- These two have a few distinct characteristics com- 
pared to the included scattering mechanisms (mainly elastic 
and spin-flip scattering): 

• They are inelastic scattering mechanisms, i.e., they 
lead to the non-conservation of spectral currents. This 
is why these should be taken into account similarly 
as the self-consistency relation, Eq. dTOb . However, 
electron-electron scattering conserves the total energy 
and charge current, whereas electron-phonon scattering 
conserves only the charge current. 



• These scattering mechanisms provide both dephas- 
ing and energy relaxation, i.e., both their Re- 
tarded/Advanced and Keldysh parts are finite. 

• Similarly to the self-consistency relation, these scat- 
tering mechanisms make the equations for the Re- 
tarded/Advanced functions depend on the distribution 
functions //, and fj . 

The self-energies for these scattering mechanisms in the pres- 
ence of superconductivity are detailed in Refill. 

Furthermore, as we concentrate only on the diffusive limit, 
we neglect effects related to different types of elastic scatter- 
ing. 

III. SUPERCURRENT SPECTRUM AND 
NONEQUILIBRIUM ELECTRON ENERGY DISTRIBUTION 
FUNCTION 

The presence of the supercurrent-induced terms js and 2F 
in Eq. ([H} leads to the finite thermoelectric effects described 
in Sec. [IV] But before engaging to their discussion, let us take 
a look at the spectral supercurrent js and how its form can be 
employed together with a nonequilibrium distribution func- 
tion to tune the supercurrent flowing in a Josephson junction, 
or alternatively, to modify the energy distribution function. 

A. Spectral supercurrent 

If a phase-coherent normal-metal wire is sandwiched be- 
tween two superconductors, Andreev reflection at each NS in- 
terface results into a formation of Andreev bound state s 27 ' 28 . 
In the case of a clean normal metal, these bound-state energies 
depend on the phase difference <p between the superconduct- 
ing contacts, the traversal time d/vp through the normal-metal 
region of length d, and the transparency z of the NS interface. 
For a junction much longer than the superconducting coher- 
ence length, the bound-state energies are 2 ^ 



= ^arcsin y t 2 cos 2 (J^j + (1 — t 2 ) sin 2 (a) + nnj 

(13) 

Here a = kpd + 5 is the dynamical phase gathered within 
traversal through the junction, 8 depending on the phase shift 
at the interface. The characteristic property of these bound 
states is that they carry an amount of supercurrent proportional 
to the phase derivative of the bound-state energy. Therefore, 
we can define a "spectral supercurrent" via 

In the clean limit js would hence contain a sequence of delta 
peaks. In the diffusive limit on which we concentrate in this 
paper, the Andreev state spectrum becomes continuous as dis- 
order makes rise to a distribution of transparencies and times 
of flights. In this case, js can be calculated by solving Eqs. (0 



FIG. 1: (a) Three-probe structure consisting of two superconduct- 
ing terminals and one normal-metal terminal. The phase difference 
<p between the superconducting terminals drives supercurrent Ig, and 
the voltage bias V in the normal terminal drives quasiparticle cur- 
rent lq P . (b) Andreev interferometer, consisting of a superconducting 
loop and two normal-metal terminals, connected by 5 normal-metal 
wires. The magnetic flux <I> threading the loop controls the supercon- 
ducting phase difference <p = 27rS>/<I>(). We take the relative lengths 
of the wires to be Lj/Lsns = liji kiji 3 an ^ assume the wires to 
have the same cross-sectional area A and conductivity a. In the nu- 
merics, we assume the wires to be quasi-one-dimensional, \/~A <C L. 
The absolute size of the system controls the characteristic Thouless 
energy scale Ej = hD/L SNS , with L S ms = ^3 +£4 +£5. 



with proper boundary conditions. Its behavior in different lim- 
its is detailed in Refsi 30 i 31 . An example of js(E) specific to 
the geometries considered in this paper is presented in Fig. [2] 

If no dc voltage between the superconductors is applied, the 
supercurrent between them is obtained from Eq. ( fTTt . 
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AEf L {E)j s {E). 



(14) 



Attaching normal-metal terminals to the wire allows one to 
tune the energy distribution function fiiE), and thereby the 



supercurrent 



30.32.33.34 



Such nonequilibrium supercurrent was 



experimentally demonstrated around the turn of the century by 
many group a 35 ' 36 i 37 ' 38 i 39 ' 4Q . One of the most interesting fea- 
tures of these experiments is the possibility to take the junc- 
tion into the 7£-state, where the ground state of the junction 
corresponds to a phase difference of % between the contacts, 
and the supercurrent for a given phase difference is reversed 
compared to the usual 0-state. 35,41 . This ft-state occurs when 
the distribution function /l weighs the negative part of the su- 
percurrent spectrum more than the positive part (c.f., Fig. O. 



B. Driving a nonequilibrium energy distribution with 
supercurrent 

Let us consider the solution to the kinetic equations ([8]) in a 
three-probe system depicted in Fig.QJa). The two supercon- 
ducting terminals are assumed to be at zero potential, whereas 
the normal-metal terminal is at potential V. For simplicity, let 
us assume the system left-right symmetric. In this case, the 
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FIG. 2: Spectrum of the supercurrent in wires 3, 4, 5 in the structure 
of Fig.[Hb), for phase difference (p = 1.6, and superconducting gap 
|A| =20E T . 



following symmetries apply inside the horizontal wire: 

js(<p) = -js(-<P) 
^{(p,x) = -£7(-(p,x) = -&{(p,-x) 
@ T ((p,x) = 9 T (-(p,x) = 9 T (<p,-x) 
@l(<P,x) = £>l(-<P,x) = @ L (<p,-x) 

In the vertical wire, we hence have j$ — £?= 0, and the kinetic 
equations for fj and fi are decoupled. Let us now try to solve 
for = fi + S/l(x) in the horizontal wire. Here ff> = 
[tznh{{E + eV)/(2k B T)) -tanh((£ - eV)/(2k B T))]/2 is the 
longitudinal distribution in the normal terminal. Using the 
fact that for < |A|, F^f = throughout the normal-metal 
system, we can find an exact solution for these energies: 



S/l(x) 







(15) 



This solution can now be substituted to Eq. (ISab . The lat- 
ter yields a second-order linear differential equation for fx, 
independent of //,. From the full numerical solution we 
can find that the proximity corrections to fj are relatively 
small compared to those in 8fi. Therefore, let us neglect 
those corrections and solve Eq. d8at in the incoherent limit 
&t = 1, ST = js = 0- in mis case we S et fr(x) = (1 ~ 
i")/r> where f T ' = PaJj is the transverse function at the 

crossing point x = 0. Here f$ = [tanh((£ + eV)/(2k B T)) + 
tanh((£ — eV))/(2ksT)]/2 is the boundary condition for fj 
in the normal reservoir, = (gvAvLsns) / {OvAvLsns + 
4-GsnsAsnsLv), and A SNS / V are the cross sections and (Jsns/v 
the normal-state conductivities of the horizontal and vertical 
wires, respectively. Substituting this solution to Eq. (TT3b fi- 
nally yields 



Sf L (x) - -f T pA 



dx 1 



Lsns Jo @l(x!) 



-Js I dx' 
10 



2x 
Lsns 



(16) 

We thus find that the supercurrent controls the antisymmetric 
part of the distribution function: for a vanishing phase gra- 
dient across the wire, 5/z, = 0. For k B T <C eV, fj defines a 
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SNS 



FIG. 3: Electron distribution function f(x,E) = i[l — fj(x,E) — 
fi(x,E)] between the two superconducting terminals in Fig.QJ. The 
bias voltage is chosen V = 30Er/e, temperature is T = iEj /kg, and 
|A| 3> Ej and <p = jr/2 are assumed. The low-energy (E ~ Ej) per- 
turbation in / arises from the L/T mixing in the proximity effect, and 
the 2V-step from the Andreev reflection, see text. 



window of energy E S [— eV, eV] in which the correction is fi- 
nite (there, ^wl, whereas /° w for |£| > |eV|). Close to 
the crossing point x = 0, Si ~ 1, and the energy dependence 
of 8/l(x) reflect directly those of 3Xx) and jg. Close to the 
NS interface jc — > ± i , tends to zero, and both of the terms 
in Eq. ( Il6l i diverge. However, their sum stays finite and the re- 
maining part is roughly proportional to the spectral supercur- 
rent jg. The full distribution function f(E,x) in the horizontal 
wire is plotted in Fig. [3] for one example value of the phase 
difference. The supercurrent-induced changes in the nonequi- 
librium distribution function were recently measured 4 —, and 
the results were in a fair agreement with the theory sketched 
above. 

The longitudinal distribution function (the energy mode) fi 
describes the response of the electron system to changes in 
the temperature 43 . In this way, the above changes in can 
be understood as supercurrent-driven modifications in the lo- 
cal temperature 44 : due to the antisymmetry of about 
the crossing point x = 0, one of the horizontal arms heats up, 
and another one cools down. Such a setup thus resembles a 
Peltier-like system. However, in this case one has to deal with 
an effective temperature T t g (for its definition, see Refs i 44,45 ), 
and it turns out that for this symmetric system the increase in 
r e ff due to the Joule heating is always larger than the changes 
due to the supercurrent. Both of these issues are settled be- 
low when considering the properties of an arbitrarily shaped 
four-terminal interferometer. 



IV. MULTI-TERMINAL THERMOELECTRIC 
COEFFICIENTS 

In this section, we apply the theory formulated in Sec. [II] 
to calculate the multiterminal transport coefficients defined in 
Eq. ([T]). Main emphasis is on the appearance of thermoelec- 



tric effects, which originate from the same mixing of the L 
and T modes that in Fig. [3] modifies the shape of the elec- 
tron distribution function. Below, we calculate all thermoelec- 
tric transport coefficients in the same example setup shown in 
Fig-Eb), a typical instance of an Andreev interferometer. The 
interference effects due to superconductivity are tuned by the 
magnetic flux <t> threading the superconducting loop, which 
adjusts the superconducting phase difference <p, and observed 
by measuring various transport properties of the wire between 
the two normal terminals. We assume here the structure to be 
left-right asymmetric, not to miss certain effects that vanish 
in completely symmetric structures. 



A. Spectral thermoelectric matrix 

Based on the above discussion, one could examine transport 
in proximity structures simply by solving the Usadel equa- 
tions numerically and evaluating the current-bias relation for 
all necessary values of temperatures and voltages at the reser- 
voirs. However, for the proximity effect, it is possible to sepa- 
rate the biases from the full non-linear response of the circuit 
by making only mild assumptions. 

First, one can note that the only part of the above equations 
that is nonlinear in the electron distribution functions / is the 
self-consistency equation (TTOb . Neglecting it is often a good 
approximation if the terminals are large compared to the rest 
of the system. Disregarding Eq. ( [Tol l, the linearity in / di- 
rectly allows one to write the charge and thermal current I' c 
and I' E entering a given reservoir i as a linear combination of 
the distribution functions f ] a {E) in all reservoirs: 20 



1 = 



I'e = 



(17a) 



dEE^ip(E)fi(E). (17b) 
PJ 



Similar decomposition has been used in the literature mostly 
for describing charge transport. 46 - 47 Below, we call the set of 
functions L'^(E) the spectral thermoelectric matrix, because 
the thermoelectric linear-response coefficients are related to it 
in a natural way: 



t'J - 

11 ~ 2k B T 

L" - 
21 ~ 2k R T 



If f E 

dEL i -l T (E)sech 2 



2k B T 



dEEL'l T (E)sech 2 



2k R T 



L'L = — — [ dEElM, (£)sech 2 ( — — 
12 2k R T J TLK ' \2knT 



(18a) 
(18b) 
(18c) 



*-^h E2 ^ E ^ h2 (^r)- (18d) 



rlj 



In principle, the functions L' J a p(E) are a generalization of the 
plain linear-response coefficients. 

The matrix element L'^g (E) can be defined explicitly as the 
a-mode current flowing in terminal i in response to a j3-mode 
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FIG. 4: Elements of the spectral thermoelectric matrix L'£„ (E) as- 
sociated with the normal terminals, i,j = 1,2, in the structure of 
Fig. |TJb). Phase difference is assumed to be (p = 1.6 and the su- 
perconducting gap |A| = 2QEj. 



FIG. 5: Elements of the spectral thermoelectric matrix L'^(E) as- 
sociated with excitations in the superconductor, i = 1,2, j = 3,4. 
Assumptions are as in Fig. [4] 



unit excitation in terminal j, at energy E: 



L'^(E)= <LS?n-r a Y J ' p (E), 



(19) 



where S?i is the surface of the z":th terminal and n the corre- 
sponding normal vector. The two-component characteristic 

potential y/j-P = , wl'^) * s assumed to satisfy the kinetic 
equations together with their boundary conditions, with the 
distribution function f' a in each terminal replaced by 8 a p 8jj. 

Examples of the energy dependence of the L'^(E) func- 
tions for the four-terminal setup in Fig. [TJb) are shown in 
Figs-SlEl The two characteristic energy scales for these co- 
efficients are, similarly as for the spectral supercurrent, the 
Thouless energy Ej — 1iD/L 2 SNS and the superconducting en- 
ergy gap A. Note that since our theory is limited to static 
situations, only L-mode (temperature) bias can be applied to 
the superconductors if they are at internal equilibrium — for 
many phenomena, the coefficients in Fig. @] are more relevant 
than those in Fig. [5] However, a nonequilibrium T-mode bias 
could be generated within the static model by inducing charge 
imbalance in the superconductors, for example by injecting 
current from additional normal-metal junctions. 

Semi-analytical expressions for the coefficients L'^ (E) can 
be found by solving Eqs. dHJ up to first order in j$ and ST. In 
systems that can be considered as a circuit of quasi- ID wires, 
this leads to a circuit theory for the distribution functions. Be- 
tween two nodes with distribution functions f l = (fj, , f[ ) and 
f 2 = (/r i fl )> one finds an expression for the spectral currents 



r/~(M- 1 -f/f 2 
is. 



(20a) 



+ f Ti(/ 2 + / 1 ) + ^(;1 + ^ 2 ), 



where t\ and £2 are Nambu spin matrices, and 



M = diag(Mr,M L ), M a = I dx^ a ( x y 

r./W)^" 1 ^"', 
Jo K ' M L M T 



y= I dx dx' sgn(x — x') 
Jo Jo 



,^ L ( x)- l ® T (x')- 1 
M L M T 



(20b) 
(20c) 
(20d) 



If node 1 (or node 2) is at a clean interface to a bulk su- 
perconductor at E < |A|, one can use the asymptotic behav- 
ior @l(x) = const, x x 2 + ^(jt 3 ), 9Xx) = j s x + 0{x 2 ) to find 
My 1 =0, t = 0, 7=±1. Using conservation of the spec- 
tral current if at the nodes and suitable boundary conditions, 
one can in this way find an approximation to (E) for any 
given circuit. The quality of this approximation is usually 
quite good — in Fig. [4] such approximations are shown with 
black lines, which almost coincide with the numerical results. 
However, the spectral equations need still to be solved to find 
out the proximity-modified diffusion constants D3l a , 3 'and 
the spectral supercurrent j$. 
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B. Symmetry relations 



As discussed in Section UTaI the normal-state thermoelec- 
tric transport coefficients are usually coupled together by On- 
sager's reciprocal relation L,'^ (B) = L J ^ a (-B) under the re- 
versal of the magnetic field. The question now is: do the 
thermoelectric coefficients induced by the proximity effect 
follow this same relation, and what else can we say about 
their symmetries. In the framework of scattering theory, it 
turns out that the Onsager reciprocity applies also in hybrid 
normal-superconducting systems. 9,48 Moreover, within the 
Usadel theory, it has been shown that the off-diagonal coeffi- 
cients L\2, L21 are always odd functions of the magnetic field 
B, whereas the diagonal coefficients Z41, Z42 are even . 15 i 18 ' 20 
Below, we review the symmetries present in the Usadel frame- 
work. 

That a form of Onsager's reciprocal relation applies for the 
Usadel model can be seen from the structure of the kinetic 
equations ([8]) and symmetries of the coefficients (O under the 
reversal of the magnetic fields B (i.e., change of sign in the 
vector potential A and the superconducting phases <p, %). The 
crucial observation is that the differential operator u in the 
kinetic equations (O, vf = 0, is related to its operator adjoint 
by22 



t _ 



= (-v) 

2|A|, 


&{-B). 



(-V) + (-V)-M (21) 



-V-J5 





Here we exploited the symmetries 2> a (—B) — $t a (B), 
3(-B) = -Mfi), 8$(-B) = &{B), and j s (-B) = -j s (B) of 
the kinetic coefficients ©• From the above relation, it follows 
that for any two-component functions <j>, p, 



Ja 



= I dYhJ, 

da 



(22) 



where the flux J — p f f (B)<j> - <j) f t {—B)p - ;sp t fi0 is what 
is left over from the integration by parts on the left-hand side. 
Especially, this flux is conserved when satisfies the kinetic 
equations for +B, and p for —B. Making now use of the 
functions applied in Eq. ( fT~9b and substituting <j) = \j/j'P(+B), 
p = y/' a (—B), the conservation of J in the volume Q. of the 
structure implies 



0= / df[p t ^-0 t ^ t p] 



d,Yh-J 



da 



(23) 



d,Yn-T a {B)(j)- dyn-r p {-B)p, 



when both i and j refer to normal terminals. In this case the 
last term in J, being proportional to jg, vanishes on the termi- 
nal surfaces and 5?j. Other terms vanish due to the bound- 
ary conditions assumed for the y/ functions. Comparing this 
result to Eq. ( fT9b . one finds for i, j referring to the normal 
terminals 



mjE,B)=V'(E-B), 



(24) 



which is a form of Onsager's reciprocal relation. 

A second class of symmetries arises from the way the co- 
efficients L'^p (E) were defined in Eq. ( fT71 ). Namely, we must 
require that 

£Z2l(£)=0, (25a) 

j 

^Lj L (E)~0 for normal terminal 2, (25b) 

; 

so that no net energy current flows to any terminal at equilib- 
rium for any temperature, and that the same applies for the 
charge current entering the normal terminals. 

The third symmetry relation is important for the thermo- 
electric effects, and is specific to the quasiclassical theory. 
Namely, if Green's function G\ is a solution to the Usadel 
equation for vector potential A and self-energy X\ [G{\, 



[V-/eAt 3 ,Gi[V-/eAf 3 ,Gi]] = [Xi[Gi],Gi 



(26) 



then, the electron-hole transformed Green's function G2 = 
— T\G\ti is a solution to the same equation for —A and self- 
energy 



X 2 [G 2 ] = -tiXi[-f 1 G 2 ti]fi. 



(27) 



For X\ [G] = —IEtt, + A[G] + ^-i^Gz^ used above, we note 

that X2 (B) — X\ (-B) — the two functionals coincide. Hence, 
the transformed Green's function describes the same physi- 
cal situation, but with an inverted magnetic field. Since elec- 
tric potentials and charge currents also change sign under this 
transformation, one finds thati2i22 



L^(E,-B) = (-l 



(28) 



This symmetry makes the off-diagonal thermoelectric coeffi- 
cients odd functions of the applied magnetic field, which is 
not in agreement with all experiments. We discuss this dis- 
crepancy in more detail in Section[VT]and in the Appendix. 



1. Charge imbalance in superconducting loops 

Below, one of the aims is to model qualitative features of 
charge imbalance in superconducting loops (see Sec. Ill CI and 
Fig. Q]) without solving the Usadel equations inside supercon- 
ductors. For this, we need some effective boundary condi- 
tions to enforce at the NS interfaces instead of the usual ter- 

Consider a supercon- 



minal assumption [case a) in Sec. Ill CI 
ducting loop with a large normal-state resistance but long in- 
elastic relaxation length [case b) in Sec. Ill Cj . Deep in the 
superconductor, we then assume that the charge current is 
carried only as supercurrent with the (BCS) spectral density 
js S(E — |A|). Due to the large resistance, we can also as- 
sume Tl<P = and f>0 = for E ^ |A|, for any solution <j) of 
the kinetic equations. Near the interface, supercurrent conver- 
sion occurs and the 5-peak in Fr<p broadens, which needs to 
be handled correctly to preserve Onsager reciprocity. Equa- 
tion ( f22b defines a flux J that is conserved in the superconduc- 
tor. By our assumptions, 7 = deep in the superconductor, for 
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E ^ |A|. The exact solution / of kinetic equations (0 thus sat- 
isfies^./ = \i/ T f T (B)f - f T f T (-B)\i/- j s (\j/ T f L + f T \j/ L ) = 
and Tif = at the NS interfaces of the loop, for any y/ that 
satisfies v{— B) \j/ = 0, regardless of boundary conditions. 

The only linear boundary condition consistent with the 
above is Yj-f = GT{\B\)fj + jsfu where Gj describes con- 
ductances related to the supercurrent conversion. For simplic- 
ity, we then assume Gj = °° at E < |A| and Gj = at E > |A|, 
which results to 



r L f = o, 
r L f = o, 



fr=0, 
r rf = jsfh , 



E<\A\ 
E>\A\ 



(29a) 
(29b) 



This acknowledges the fact that for E < |A| the kinetic equa- 
tions imply a vanishing fj beyond the current conversion re- 
gion, and that in a BCS superconductor fj does not relax 
at E > |A| if there is no inelastic scattering. 22 Employing 
Eq. d29b is analogous to requiring that the "non-equilibrium" 
parts of the spectral currents vanish; the remaining part jsfi 
is what at equilibrium gives rise to the supercurrent. 

Note that Eq. d29l ) is not exact: we at least neglect the resis- 
tance in the supercurrent conversion region discussed for ex- 
ample in Refs i 22 i 49 . Note also that when treating a supercon- 
ducting loop as two boundary conditions, charge conservation 
must be ensured by adjusting all potentials relative to that of 
the superconductor. Nonetheless, we expect that Eq. d29l cap- 
tures some of the relevant physics in the problem. Below, we 
use it to illustrate how charge imbalance could change observ- 
able quantities. 



C. Conductance 

How the proximity effect changes the conductance has 
been studied in detail, both experimentall y 48 ! 50 ' 51 ! 52 ' 53 ! 54 ! 55 
and theoretically 46,56 ' 5758 . For a review, see for example 
Ref^. 

The modification to conductance can conveniently be de- 
scribed with the Usadel equations. Once L 1 ^ (E) is known — 
usually the zeroth order in j$ and 3? is accurate enough — 
calculating various conductances can be done: one can di- 
rectly evaluate the corresponding conductance matrix Ljj and 
thermoelectric coefficients L l ( 2 from Eqs. cfT~5T > and write 

d4=£4d^+£4d7>/r 

i i 



y dii_ 



(30) 



d<pj 



The second sum is finite if the heating of the terminals is sig- 
nificant, but should still give only a small contribution as the 
thermoelectric coupling is small, as can be seen in Fig. [4] 
The last term arises if conductances are evaluated in struc- 
tures where the phases <pj in the superconducting terminals 
may vary. However, for i referring to a normal terminal, 
F C ({V} — 0,{<p}) = independent of the phases {<p}. This 
implies that the last term vanishes for conductances around 




FIG. 6: Linear-response electrical resistance R between terminals 1 
and 2 of the structure in Fig. \V[b), as a function of the phase dif- 
ference <p and the temperature T . The resistances are normalized to 
the normal-state resistance — R\ + R2 +R5. The curves corre- 
spond to different models of the superconducting loop discussed in 
Sections lilCllIVB II - a) (solid) and b) (dashed). The current flows 
via the superconducting loop as supercurrent, reducing the resistance 
from the normal-state value also at T = 0. Temperature dependence 
of the energy gap A is neglected, and we assume A = 2QEj- 



{V} = 0, the potential of the superconductors, but it may be 
finite when calculating differential conductances. Note also 
that when modeling superconducting loops using only bound- 
ary conditions at the NS interfaces, current conservation needs 
to be ensured by adjusting all potentials relative to that of the 
superconducting condensate. 

Typical behavior of conductance in an Andreev interfer- 
ometer is illustrated in Fig. [6] The proximity effect adds an 
enhancement that oscillates with the superconducting phase 
difference <p and has a re-entrant dependence on the tempera- 
ture T. The figure also shows how charge transport via quasi- 
particles (E > |A|) in the superconducting loop may change 
the conductance at high temperatures. The two curves corre- 
spond to the terminal a) and long-loop b) limits discussed in 
Sections III CI IIVB II For the former, the loop contributes to 
electric conduction at energies E > |A|, for the latter it does 
not. 



D. Thermal conductance 

As for the electrical conductance, the proximity of super- 
conductors modifies also the thermal conductance i 60 ' 6 1 ' 62 ' 63 
This was studied on the basis of the quasiclassical Usadel the- 
ory in Refs^iH. 

For a given setup, calculation of the thermal conductance 
from L'^p (E) proceeds as for the electrical conductance. Typ- 
ical predicted features are ^-periodic suppression of thermal 
conductance at low temperature ksT < Ej due to modified 
density of states and thermal diffusion coefficient D&i,, and 
inhibition of sub-gap thermal transport into the superconduc- 
tors due to Andreev reflection. These are illustrated in Fig. [7] 
for the example setup, together with two models for the above- 
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FIG. 7: As Fig. [6] but the thermal resistance R t ^ is shown. It 
is normalized to the normal-state Wiedemann-Franz value R±.n — 
3e 2 R^/(7t 2 kBT). The two curves correspond to same models for the 
superconducting loop as in Fig.[6] a) terminal (solid) and b) long loop 
(dashed). In the former, at kgT ~ A, part of the thermal current flows 
through the loop as quasiparticle excitations, reducing the thermal 
resistance. Note that the scale for R,h is the same in both figures. 
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FIG. 8: Linear-response thermopower in the structure of Fig.[TJb), 
as a function of the phase-difference <p and the temperature T. 
Solid line: no charge imbalance in superconducting loop [case a) 
in Sec. Ill CI . Dashed line: no inelastic relaxation in the long super- 
conducting loop [case b) in Sec. Ill CI . Dotted line: Approximation 
{33), neglecting contributions from S. If the terms proportional to 2? 
are taken into account, the result coincides with the solid line. Other 
assumptions are as in Fig. [6) 



gap quasiparticle transport in the superconducting loop. 



E. Thermopower 

Thermopower S is proportional to the upper right coeffi- 
cient Z42 of the thermoelectric matrix. The superconducting 
proximity effect on S has recently been studied experimen- 
tally, see Refs, 64 i 65 i 66 i 67 ' 68 i 6 V° Theoretically, predictions for 
the thermopower in hybrid normal-superconductor structures 
have been calculated starting from the scattering theory infill, 
and via the Usadel theory discussed here i 15 i 16 i 17 i 18 i 19 We dis- 
cuss the comparison between theory and the experiment in 
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FIG. 9: As Fig. [8] but showing the NN thermopower Saw- The dotted 
line includes the terms proportional to ^in Eq. ( 1331 : the other terms 
in Eq. i33\ vanish. 



Section[Vl] and consider here only the theoretical model. 

For a two-probe structure, the thermopower is usually de- 
fined as the induced voltage divided by the temperature dif- 
ference when no charge current flows, S = L Q , but the 
additional terminals in the 4-probe structure in Fig.QJb) allow 
for defining two distinct thermopower-type quantities, 



_ d(V!+V 2 ) 
2>NS = 



2d(7j - T 2 ) 



_d(Vi-V 2 ) 
Saw = 



lr 1 — Ir 



d(7i - T 2 ) 



Ir 1 —Ir 



(31a) 



Both of these can be calculated from L^„ 



(£): 



-[12k_i,[12] / 1 




(32a) 
(32b) 

(32c) 



Typical results are shown in Figs. [8] and [9] The oscillations 
in <p are always antisymmetric due to the symmetry rela- 
tion (128k and the temperature dependence shows the reen- 
trant behavior on the energy scale of Ej characteristic of the 
superconducting proximity effect. One can also note that 
the magnitude of the effect is significantly larger than what 
is expected from the normal-state thermoelectric effects at 
sub-Kelvin temperatures, which typically are of the order of 
S« 10~ 4 ... 10~ 3 xk B /e. 

Making use of expression ( f20l > and neglecting the energy- 
dependence of £>t and @l one can aiso derive approximations 
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to the charge current at constant temperature, II = . In our 
example four-probe structure in Fig. |TJb), two Peltier coeffi- 
cients can be defined, 



FIG. 10: As Fig. [8] but showing the Peltier coefficient Uns- The dot- 
ted line is obtained from approximation J33t including the 2T~ terms. 
The Kelvin relation TL = TS can be seen by comparing to Fig. [8] 



such as 17 



Snn ■ 



(R 3 - 



R 4 )R 2 5 



dl s 



cq 



2(R l +R 2 + R 5 )R SNS dT 
RsNs(bi+b 2 ) + {R 3 +R 4 )b 5 



+ 



4R 3 R 4 R 5 -f 



\-R 2 + Rs)Rsns 
R 2 5 {R 3 +R 4 )dI s 



eq 



4(R 1 +R 2 +R 5 )Rsns dT 
R S Ns{bi-b 2 ) + (R 3 -R 4 )b 5 



(33a) 



(33b) 



2(R l+ R 2 +R 5 )R, 



SNS 



where R SNS = R 3 +R 4 +R 5 , |A| > E T , and 



dEE 



2eksT- 



■ sech" 



( E 



\2k B T J Lj Jo 



R, 



dx^x) (34) 



are averages of the coefficient 3? in different wires. The ap- 
proximation (l33l is compared to the numerical solution in 
Figs. [8] and [9] It turns out that a large part of the ther- 
mopower is related to the equilibrium supercurrent Is,eq> 15,17 
Note also that the contribution from Is to Snn 1S strongly de- 
pendent on the asymmetry in the structure and vanishes for 
a left-right symmetric setup, as does the contribution from 
ffi 17 ' 18 However, the contribution from energies E > |A|, 
which is neglected here, behaves differently in this respect, 
see Refs»i2. 



F. Peltier effect 

The second off-diagonal thermoelectric coefficient L 2 \ has 
not yet been measured in the presence of the proximity ef- 
fect, although related experiments far from equilibrium have 
been made. 42 Theoretical predictions for modifications due to 
the proximity effect have been calculated from the scattering 
theory 9 and from the Usadel theory^ -. 

A finite L 2 \ coefficient induces a Peltier effect, energy cur- 
rent driven by charge current. The Peltier coefficient IT is in 
general defined as the ratio of the heat current Iq=Ie — V-h 



IT 



NS 



d4 

dl c 



n 



NN 



d}± 

dl c 



I}=-I}=lc 



(35) 



corresponding to two different current configurations. These 
are directly related to the linear-response L-coefficients by 



(36a) 
(36b) 



in a similar way as in Eq. ( f32t . However, note that IFvs can be 
defined only when there is a grounded extra contact in the su- 
perconducting loop [c.f. Fig. [Tib)] through which the injected 
current I c can flow. 

As discussed above, the matrix element L 2 \ is usually cou- 
pled to the element L\ 2 via Onsager's reciprocal relation. This 
leads to Kelvin relations between the Peltier coefficients and 
the thermopower 



n#s — tsns , 



FIaw — TSnn , 



(37) 



which are easily seen by transposing equations d36l ) and 
comparing to Eqs. d32| l. These relations are not broken by 
the superconducting proximity effect, which implies that the 
proximity-induced Peltier coefficient inherits the magnitude, 
phase oscillations and the temperature dependence of the ther- 
mopower. Numerically calculated linear-response Peltier co- 
efficient in the example structure is illustrated in Fig.QJJ 

The Peltier coefficient is sufficiently large so that it could 
be detected simply by observing how the effect changes the 
temperature of one of the terminals in Fig. |TJb). For a 
typical Thouless energy Ej/ks = 200 mK, the coefficient in 
Fig. [10] achieves a magnitude of IT ~ 1.5 jxW at temperature 
T = 400 mK. A simple heat balance estimate, assuming that 
the terminal 1 is thermally isolated apart from the electronic 
heat conduction through wire 1, 



Iq = —G t hAT- 



2U NS I c + eI 2 /G = 0, 



(38) 



then 



yields 

-(3/^ 2 )( e 2 n^ 



a maximum cooling AT w 
jk\T) ~ 0.2 mK. However, the oscilla- 



tion amplitude is proportional to I c and can be larger than this 
maximum cooling effect: variation of the order of millikelvin 
at least should be possible. 20 Temperature changes of this 
order have already been successfully resolved in mesoscopic 
structures 71 , so that in a suitably optimized setup, it might 
also be possible to detect this proximity-Peltier effect. 



V. DEPENDENCE ON EXTERNAL FLUX 

A magnetic field applied to a normal-metal-superconductor 
heterostructure causes persistent currents to flow in the struc- 
ture and induces some dephasing. The currents also screen 
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the applied magnetic field, which can usually be taken into 
account by assigning self-inductances to all loops in the struc- 
ture. Both effects can be included in the present theory, and 
we discuss the latter briefly below. 

If considering the Andreev interferometer in Fig. 02b), 
screening is mostly taken into account in the / c (<p) relation 
of the weak link. The inductance L of the loop only modi- 
fies the <p{<b x ) relation between the induced phase difference 
(p and the external magnetic field <t> x to 26 ' 70 



4> 

q>-2n—=U c (q>). 



(39) 



One should note that although a modified <p(& x ) relation 
should change the shape of the oscillation of various quantities 
as functions of 4> v , e.g. thermopower in Fig. [8] the symmetry 
properties in Sec lIVBl remain unchanged. However, if there 
is hysteresis and multiple flux states are possible for the same 
values of control parameters, the situation is slightly more 
complicated: for a given solution of ( 1391 with external flux 
<& x , there exists a solution with — for which Eqs. (124128b 

apply- 
There is a further effect of the magnetic field neglected in 
this work: Zeeman effect, which leads effectively to an ex- 
change field inside the wires (for an example of such an effect, 
see Ref. 72 ). However, unless special care is taken, this effect 
plays typically a much smaller role than the dephasing effect 
of the field. 



VI. DISCUSSION 

In this article, we have systematically discussed the predic- 
tions of the quasiclassical diffusive-limit theory on the ther- 
moelectric response of normal-metal samples under the influ- 
ence of the proximity effect. The latter yields corrections to 
the fairly general relations in Eq. ©. These corrections de- 
pend in general on energy (i.e., on temperature or voltage) and 
on the phase difference between superconducting contacts. At 
least in most typical cases, one of the general relations, the 
Onsager relation (and thereby also the Kelvin relation) holds 
also in the presence of the proximity effect. Furthermore, 
the approximations made in the quasiclassical theory imply 
that the diagonal coefficients of the thermoelectric matrix are 
generally symmetric and the off-diagonal ones antisymmetric 
with respect to an external magnetic flux. 

Our results for the proximity correction of the conductance 
agrees with the previous quasiclassical treatment s - 46 ' 56 ' 57 ' 58 . 
However, as far as we know, the charge imbalance effect has 
not been previously addressed. The thermal conductance cal- 
culated here is in line with the results in Ref. 67 , but in contrast 
to it, we do not make any approximations to the kinetic equa- 
tions. 

The quasiclassical prediction on the thermopower has been 
detailed in different situations in Refs. 1 5 ' 16 ' 17 ' 18 ' 19 . Our the- 
ory is in line with these predictions. The mechanism for the 
finite thermopower is analogous to the generation of charge 
imbalance in bulk superconductors in the presence of coexist- 
ing supercurrent and temperature gradien t 13 ' 73 ' 74 ' 75 . 



To our knowledge, the only quasiclassical treatment of the 
Peltier effect and the resulting temperature modification prior 
to this paper is our Ref. 20 . Beyond the quasiclassical approx- 
imation, these effects have been discussed using the scatter- 
ing theory and numerical simulations of the Bogoliubov - de 
Gennes equation on a tight-binding lattice 9 -^. In that work, 
the symmetry of the flux dependence for the off-diagonal co- 
efficients was mostly dependent on the geometry and disor- 
der of the considered system, and not fixed as in our work. 
However, the small size of the simulated structures makes a 
quantitative comparison for example to the present work dif- 
ficult: in Ref* 9 - even the normal-state thermoelectric effects 
were large, and it is difficult to distinguish those contributions 
from the proximity effect that remain large in experimentally 
relevant structures from those that rely on significant electron- 
hole asymmetry. 

On the experimental side, a qualitative agreement to most 
of the features presented here has been found. The resistance 
correction in an Andreev interferometer has been found to os- 
cillate with a magnetic flux through the loop 48,50 - 5 1 i 52 i 53 t 54 i 55 , 
with the scale given by the flux quantum. Moreover, the reen- 
trance effect illustrated in Fig. [6] has been measured in dif- 
ferent sample s 48 ' 53 ' 55 . However, to our knowledge there is 
no successful quantitative fit between the quasiclassical pre- 
dictions and the experimentally measured temperature depen- 
dence of the resistance — see an example of such a com- 
parison in Ref. 55 . The reason for this may be the neglect 
of the generally temperature-dependent inelastic scattering ef- 
fects (see Sec. Ill Dl ) in the theory. 

The thermopower in the presence of the proximity effect 
has been measured by two groups, one in the Northwestern 
University, USA 64,65,66,67 , and another in the Royal Holloway 
University of London 68 - 69,70,76 . Again, most of the qualitative 
features agree with the quasiclassical theory. The measured 
thermopower oscillates with the flux and is at least two or- 
ders of magnitude larger than the normal-state thermopower, 
and in line with the predictions from the quasiclassical theory. 
The first attempt for a quantitative fit 70 of the temperature de- 
pendent thermopower between the theory and the experiments 
was unsuccessful. We believe that the major reasons for this 
were the too complicated geometry of the measurements for 
this purpose and the neglect of the inelastic scattering effects. 

The major qualitative disagreement between the theory 
and the measurement is in the symmetry of the ther- 
mopower oscillations with the flux: in most measurements, 
the oscillations were antisymmetric and in line with the 
theory 64,65,66,67 ,6 8 i 69 ' 70 , in some measurements they were 
symmetric 66,67 . The authors of Ref. 67 suggested that this 
symmetry depends on the geometry of the sample: in sam- 
ples where the supercurrent flows along with the temperature 
gradient, the oscillations are antisymmetric whereas in other 
types of samples they are symmetric. Such a conclusion can- 
not be made based on the quasiclassical theory. 

We also note that in bulk superconductors, the magnitude 
of the thermoelectric effects has been long under debate 77 — 
there the experiments have shown larger thermoelectric effects 
than those predicted by the theory. 

The only published measurement on the thermal resistance 
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/?th of an Andreev interferometer known to us» showed an 
oscillating R t h, but the correction from the proximity effect 
was larger than that predicted by the theory. We are not aware 
of any measurements of the Peltier effect. 

Quasiclassical theory, based on the combination of the BCS 
model and the quasiclassical approximation, has been suc- 
cessful in providing a quantitative explanation to a broad 
range of superconducting phenomena. Here we have pointed 
out one qualitative aspect (flux symmetry of the thermoelec- 
tric effects) which is yet to be explained. Clearly, the full un- 
derstanding of the nonequilibrium electron transport phenom- 
ena in superconducting proximity samples will still require 
both further experimental and theory work. 

APPENDIX A: POSSIBLE REASONS FOR THE 
SYMMETRIC THERMOPOWER OSCILLATIONS 

In the diffusive limit, the antisymmetric flux dependence of 
the proximity-induced thermopower results from the special 
symmetry of the self-energies: all the typically relevant self- 
energies satisfy Eq. (|27l i in the presence of a magnetic field 
B with X2(B) = X\ (-B). Outside the diffusive limit, one has 
to employ the Eilenberger equation 78 describing the Keldysh 



Green's function g(p,r,E ,B). Here p is the direction of the 
electron momentum and r is the center-of-mass coordinate. In 
this case, the property of the self-energies x[g] leading to the 
antisymmetric thermopower oscillations is 

x[g(p,r,E,B)} = -Til[-tig(-/3,F,£,-fi)ti]Ti. (Al) 

This symmetry is satisfied for the most relevant self-energies, 
including those for the elastic or spin-flip scattering in the 
Born approximation, and that related to the superconducting 
order parameter. We note that in Ref.— , it was shown that a 
dilute concentration of impurities away from the Born limit 
leads to large thermoelectric effects in unconventional super- 
conductors. 

Beyond the quasiclassical approximation, other possible 
reasons for the symmetric thermopower oscillations may be 
largely enhanced electron-hole asymmetry effects (however, 
these were shown in Refs. 80 ' 81 to be small for a fairly generic 
setup) or quantum interference contributions 82 . Further stud- 
ies on these effects are therefore required. 

This research was supported by the Finnish Cultural Foun- 
dation and the Academy of Finland. We thank N. Birge, M. 
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